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Predictions for the hyperfine splitting of 2S bottomonia are compared with the two recent ex- 
perimental candidates for the jjt(2S). The smaller splitting of the Belle state is consistent with 
unquenched lattice QCD computations, many potential models, and a model-independent relation 
which works well for charmonia. The larger splitting for the state extracted from CLEO data is 
inconsistent with most predictions. 
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I. INTRODUCTION 

Within a short space of time, two groups have indepen- 
dently claimed the discovery of the r)h(2S). In an analysis 
of CLEO data, Dobbs et al. [1] observe an enhancement 
in the Ml decay T(2S) —> rfo(2S)j, with mass and hy- 
perfine splitting AM 2S = M T(2S ) - M Vb{2 s)- 

M Vb (2S) = 9974.6 ±2.3 ±2.1 MeV, (1) 
AM 2S =48.7 ±2.3 ±2.1 MeV. (2) 

The Belle collaboration [2] reports a state in the El decay 
/ib(2P) — > ?7b(2S)7, with mass and hyperfine splitting: 



M, 



%(2S) 



9999.0 ± 3.5±f J MeV, 



AM 2S = 24.3; 



■4.o 



MeV. 



(3) 
(4) 



The striking disagreement in the masses has not yet been 
addressed in the literature. The aim of this paper is 
to compare the conflicting results with predictions from 
a range of theoretical approaches: lattice QCD (Sec- 
tion II), a model-independent mass relation (Section III), 
potential models (Section IV), and the unquenched quark 
model (Section V). 

The data which forms the basis of the discussion is 
collected in Table I, a compilation of experimental val- 
ues and theoretical predictions for both the IS and 2S 
hyperfine splittings of bottomonia. It is useful to look 
at the IS and 2S splittings together, because the two 
are not independent: several of the lattice QCD pre- 
dictions for the 2S splitting, and the prediction of the 
model-independent mass relation, are normalised against 
the IS splitting, and in potential models the IS and 2S 
splittings are controlled by the same parameters of the 
spin-spin interaction. 

The role of the IS splitting in the interpretation of 
predicted 2S splittings is particularly important because 
its experimental value is still in flux. Various mea- 
surements and averages for the IS splitting are shown 
in Table I. The 2012 PDG value ("PDG12") is ob- 
tained from an average of three experimental results in 
T(3S) -)• 7? b (lS)7 [3, 4] and T(2S) -» %(1S)7 [5]. More 
recently Belle [6] observed /i&(1P) — > f?b(!S)7 and their 
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?7fc(lS) mass corresponds to a hyperfine splitting with cen- 
tral value 10.0 MeV lower than PDG12 . In a mini-review 
in the 2012 PDG [7, 8], a new average value for the IS 
splitting ("PDG12*" in Table I) is computed taking ac- 
count of the Belle measurement, almost 5 MeV smaller 
than the previous average. 

Dobbs et al. and Belle also observe the 77b(lS) in 
the analogous radiative decays in which they claim the 
?7fc(2S); their measured values for the IS hyperfine split- 
ting are also shown in Table I. In an analysis of T(1S) — >• 
7yb(lS)7 Dobbs et al. obtain a IS splitting which is consis- 
tent with PDG12 and PDG12*. In h b (lP, 2P) -» %(1S)7 
Belle obtains a IS splitting which is smaller than both 
PDG12 and PDG12*, but consistent with the earlier 
Belle [6] measurement. 

The important point is that each of the experimen- 
tal measurements which post-date the computation of 
PDG12 have a smaller central value for the IS splitting. 
This is particularly relevant to the analysis of the lattice 
QCD predictions, which is the first topic of discussion. 



II. LATTICE QCD 

Various predictions for the IS and 2S splittings in non- 
relativistic lattice QCD are collected in Table I. Apart 
from those marked (o), all of the results are from un- 
quenched lattice QCD. 

For the IS splitting, Gray et al. [9] (pre-dating the 
discovery of the r/i, (IS)) and Meinel [10] are consistent 
with all of the measured values and the PDG averages; 
Dowdall et al. [11] and the quenched result of Burch and 
Ehmann [14] are consistent with the PDG averages and 
Dobbs et al, while Lewis and Woloshyn [12, 13] are con- 
sistent with Belle. 

The predictions for the 2S splitting discriminate more 
strongly. The result of Gray et al. [9], with its some- 
what larger errors, is consistent with both Dobbs et al. 
and Belle. Otherwise, none of the predictions are con- 
sistent with Dobbs et al. The predictions of Meinel [10] 
(one normalised to the IS splitting, the other to the IP 
splitting), Lewis and Woloshyn [12, 13], and Burch and 
Ehmann [14] are all consistent with Belle. The prediction 
of Dowdall et al. [11] falls between the measured values of 
Belle and Dobbs et al. and is not consistent with either. 

Two of the predictions for the 2S splitting, marked in 
the table (t), are normalised to the PDG12 value of the 
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Lattice QCD 
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TABLE I. Experimental data and theoretical predictions for 
the IS and 2S hyperfme splittings (in MeV) of bottomonia. 
t Normalised to the experimental IP tensor splitting, 
j Normalised to the PDG12 IS splitting, 
o Quenched lattice QCD. 



IS splitting. As discussed previously, recent experiments 
[1, 2, 6] all measure a smaller value for the IS splitting, 
which would imply that these lattice predictions are over- 
estimates. This is particularly interesting in the case of 
Dowdall et al. [11], the only lattice prediction which is 
inconsistent with Belle. 

The quantity measured on the lattice is the ratio of 2S 
to IS splittings. Dowdall et al. obtain 



2S 



0.499 ±0.042, 



(5) 



AM 1S 

while Meinel obtains (with errors added in quadrature) 
AM 2S 



AM- 



0.403 ±0.059. 



is 



(6) 



By normalising against the updated PDG12* IS split- 
ting, rather than PDG12, the prediction of Dowdall 
et al. is brought into agreement with Belle, AM 2 s = 
32.2 ± 4.2 MeV, and the agreement between Meinel and 
Belle is improved, AM 2S = 25.9 ± 5.0 MeV. 

An alternative interpretation of the results of Dow- 
dall et al. and Meinel is to normalise against their com- 
puted values for the IS splitting. The resulting predic- 
tions are consistent with Belle, AM 2 s = 34.9 ± 7.4 MeV 
and AM 2S = 24.3 ± 6.7 MeV respectively. 

It is also interesting to normalise Dowdall et al. and 
Meinel against the most accurate measurement of the IS 
splitting from a single experiment, that of Belle [2]. The 
predicted 2S splittings are consistent with the measured 
value at Belle despite smaller errors, AM 2 s = 28.9 ± 
3.6 MeV and AM 2S = 23.3 ± 4.3 MeV respectively. 

(Notice that the adjustment to the Meinel value brings 
it into better agreement with the other value of Meinel 
quoted in Table I, which is normalised to the IP tensor 
splitting. This implies that the ratio of IS hyperfme to 
IP tensor splittings, which Meinel also computes directly 
on the lattice, is in better agreement with experiment us- 
ing the IS of Belle than with PDG12, which is interesting 
in its own right. Intriguingly, the computed masses of B c 
and B s mesons in ref. [51] are in better agreement with 
experiment when normalised to the Belle ^(lS) mass 
rather than that of the PDG, although it is not a statis- 
tically significant effect.) 

A more direct interpretation of the above lattice re- 
sults is to compare the computed ratio of splittings di- 
rectly with experiment, rather than extracting a predic- 
tion for the 2S splitting which depends upon (and is sub- 
ject to the error in) the IS splitting against which it is 
normalised. In this context, two other lattice results are 
worth mentioning. Gray et al. [9] measure the ratio 



AM- 



2S 



AM 1£ 



0.5 ±0.3, 



(7) 



while the quenched computation of Davies et al. [52] has 
the ratio 



AM 2S 
AMis 



= 0.49 ±0.09. 



(8) 
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(The predicted splittings of the latter do not appear in 
Table II because they are presented in lattice units.) 
From the Belle data, one obtains the ratio 



AM 2S 
AM 1S 



= 0.42 ±0.09, 



(9) 



in agreement with all of the above. On the other hand 
the ratio of Dobbs et al. 



AM 2S 
AM- 



0.73 ±0.09 



is 



(10) 



is consistent only with Gray et al. [9]. 

Figure 1 shows a plot of the IS and 2S splittings of var- 
ious theoretical approaches compared with experimental 
data. The measures splittings of Dobbs et al. and Belle 
are shown by error bars, and the lattice predictions shown 
by the shaded areas. (The lines and circles show the pre- 
dictions of the model-independent relation and poten- 
tial models, discussed in the subsequent sections.) The 
skewed shapes are the results of Meinel and Dowdall et al. 
using their measured IS splittings and the measured ra- 
tios AM2s/AMis- The remaining two data points are 
the other result of Meinel (which is normalised to the IP 
tensor splitting) and that of Lewis and Woloshyn. For 
the sake of clarity, the results of refs [9, 14] with larger 
errors are not plotted. 



III. MODEL-INDEPENDENT RELATION 

In this section the experimental 2S splittings are com- 
pared to the predictions of a mass relation which is com- 
mon to the simplest non-relativistic potential models. 
The hyperfine splitting is controlled by the spin-spin term 
in the potential, which in tree-level perturbative QCD is 
proportional to a delta function in the quark separation, 



q • Sq, 



(11) 



where m is the quark mass and a s controls the strength of 
the interaction. In perturbation theory the correspond- 
ing hyperfine splitting AM nS = M n 3 Sl - M„i So is pro- 
portional to the square of the radial wavefunction R n s at 
the origin: 



AM nS 



9m 2 



|i?™s(0)| 2 . 



(12) 



The e + e~ width of an n 3 Si meson is also proportional to 
the square of the wavefunction at the origin, according 
to the Van Royen-Weisskopf formula [53] , 



r fi + f 



^ 3 Sl 



4e 2 a 2 



16a f 



(13) 



where e q is the electric charge of the quarks in units 
of the fundamental charge, and the a s term is the first 
QCD correction [54]. For M n one normally uses twice the 



quark mass, or the mass of the n 3 Si state; as discussed 
by Voloshin [55] the two prescriptions are equally valid 
in the non-relativistic limit. 

The model-independent wavefunction at the origin 
cancels in the ratio of hyperfine splitting to e + e~ width: 



AM 



nS 



i 3 Si 



9e 2 q a 2 



Mn 

2m 



1 — 16a s /37r 



(14) 



With the prescription M n = 2m, the experimental data 
on the IS levels of charmonium and bottomonium satisfy 
the above with reasonable values of a s = 0.29 and a s — 
0.22 respectively. Taking ratios of the above equation for 
IS and 2S levels yields a relation (which depends upon 
the prescription for M n ) among the hyperfine splittings 
and e + e _ widths of the IS and 2S levels. With M n = 2m, 



AM 



2S 



Pp+ f 



>2 3 Si 



AM 1S 

and with M n — M„3 Sl , 

AM 2S r e+e 



>l 3 Sl 



AM, 



IS 



>l 3 Si 



M 



2 3 Si 



M- 



l 3 Si 



(15) 



(16) 



Before applying the relations to bottomonia, it is inter- 
esting to test their validity for charmonia. The difference 
between the two prescriptions is quite significant, with 
(M^/Mj/^,) 2 = 1.416. The experimental data [7] gives 
the ratios 



AM 2S 
AMis 



0.407 ±0.015, 



= 0.423 ±0.018, 



e+e —¥j/ip 



(17) 
(18) 



which are consistent with the relation taking the pre- 
scription M n = 2m. 

For bottomonia, the difference between the two pre- 
scriptions is less significant, because (Mf /Mr) 2 = 
1.123. This is another manifestation of the familiar sit- 
uation that non-relativistic predictions are more power- 
ful for bottomonia than charmonia. The ratio of e + e~ 
widths [7] is 



= 0.457 ±0.014. 



(19) 



The ratio (9) of hyperfine splittings at Belle therefore 
satisfies the model-independent relation using either pre- 
scription for M n . On the other hand the ratio (10) of 
Dobbs et al. does not satisfy the relation using either pre- 
scription. The prediction of the relation (with M n = 2m) 
is shown by the broken lines in Figure 1. 

Notice also that the ratio of splittings predicted by the 
model-independent relation is consistent with all of the 
ratios (5), (6), (7) and (8) computed on the lattice. From 
Figure 1, it is also clearly consistent with the other lattice 
computations. 
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FIG. 1. Experimental data and theoretical predictions for the hyperfine splittings of IS and 2S bottomonia. The data points 
of Dobbs et al. and Belle are shown by the error bars. The predictions of lattice QCD (Section II) are shown by the shaded 
areas, the model-independent mass relation (Section III) by the broken lines, and potential models (Section IV) by circles. 



IV. POTENTIAL MODELS AND RELATED 

In this section the observed splittings at of Dobbs et al. 
and Belle are compared to the predictions of potential 
models and related approaches. Table I collects a range 
of such predictions, including ordinary non-relativistic 
quark potential models with various functional forms 
(Cornell, logarithmic, this power law, Buchmullcr and 
Tye, Richardson's, screened); scalar, vector and mixed 
confining potentials; relativistic, relativised and semi- 
relativistic models; perturbative QCD and renormalisa- 
tion group-improved potentials; and one model in which 
the potential is extracted from lattice QCD. The physics 
of the various approaches will not be discussed in detail 
here: the aim of this section is to draw more general 
conclusions as to the feasibility of fitting the disparate 
experimental results to potential model predictions. 

Some of the approaches quote several predictions for 
the IS and 2S splittings, corresponding to different pa- 
rameter sets. In such cases, if a particular parameter set 
leads to a worse fit to all of the experimental data (PDG, 
Dobbs et al. and Belle) than another parameter set, it is 
not included in the table. Those parameter sets which 
cannot be discriminated in this way are all included in 
the table. 

The data in Table I is plotted in Figure 1. There is 
considerable spread in the predictions of both the IS and 
2S splittings. There is one data point which is consistent 
with the 2S splitting of Dobbs et al. but the correspond- 
ing prediction for the IS splitting is much larger than any 



of the experimental values. These splittings (the third 
pair of values quoted in the table for Bali et al. [17]) 
are obtained by solving the Schrodingcr equation with 
a potential extracted from quenched lattice QCD, and 
incorporating an estimate of the correction due to the 
quenched approximation. The same approach applied to 
charmonia predicts IS and 2S hyperfine splittings which 
are much larger than experiment. 

All of the other models predict the 2S splitting smaller 
than that measured by Dobbs et al, in most cases quite 
considerably smaller. Those which lie closest to 2S split- 
ting of Dobbs et al. all predict the IS splitting much 
larger than experiment. 

On the other hand many of the predictions [17-21, 25, 
27, 29, 32-35, 39-42, 45, 47, 49, 50, 56] are consistent with 
the 2S splitting of Belle, among which the IS splittings of 
several [18, 19, 39-41, 47] are also consistent with Belle, 
two [29, 35] are consistent with Dobbs et al. and PDG12*, 
and the rest [17, 20, 21, 25, 27, 32-34, 42, 45, 49, 50, 56] 
are smaller than any of the measured values. 

The data in Figure 1 shows a positive correlation be- 
tween the IS and 2S splittings, as expected since both are 
controlled by the same parameter (a s ) in the potential. 
It is plausible that models for which the ratio (but not 
the magnitude) of 2S and IS splittings is in agreement 
with experiment can be brought into agreement also in 
magnitude by suitably adjusting a s . Of course a s ap- 
pears in other terms in the potential, so whether or not 
such an adjustment improves or spoils the overall fit to 
the spectrum will vary from model to model. Neverthe- 
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less it is noteworthy that the ratio of 2S to IS splittings 
of Dobbs et al. is not consistent with any of the models 
considered, as can be seen in Figure 1: all of the data 
points would lie below the line of minimal slope between 
the origin and the error bars of Dobbs et al. On the other 
hand, the measured ratio of Belle is consistent with many 
models [15, 16, 18-24, 26-30, 33-35, 37-43, 45, 47, 48]. 

Finally it is interesting to consider separately the pre- 
dictions of those models which came before and after the 
discovery of the ^(lS). All of the models which post- 
date the discovery of the %(1S) [15, 16, 28, 37, 46, 48] 
are fit to the earlier IS splittings measured at BaBar [3, 5] 
and CLEO [4], or to the PDG averages at the time, and 
do not take account of the more recent, smaller values of 
Belle [2, 6] and Dobbs et al. [1] shown in Table I. Each of 
these models also predicts a 2S splitting somewhat larger 
than that observed at Belle. Presumably an updated 
fit in these models, which takes account of the smaller 
IS splittings of more recent experiments, will decrease 
the predicted 2S splittings, possibly into agreement with 
Belle. 

Among the predictions for the IS splitting from models 
which which pre-date the discovery of the ^(lS) there 
are, in decreasing order of splitting: several which are 
larger than any of the measured values [17, 21, 23, 24, 
26, 36, 38, 43]; one [30] consistent with PDG12; two 
[29, 35] consistent with the smaller PDG12* and Dobbs 
et al. values; several [18, 19, 22, 39-41, 44, 47] consis- 
tent with the even smaller value of Belle; and quite a few 
[17, 20, 21, 25, 27, 31-34, 42, 45, 47, 49, 50, 56] which are 
smaller than any of the measured values. The implication 
is that models which are not constrained to fit the r?b(lS) 
mass typically predict a smaller IS splitting, which is in- 
teresting in light of the smaller measured values at recent 
experiments. 



V. THE UNQUENCHED QUARK MODEL 

The predictions of the model-independent relation and 
quark potential models discussed in Sections III and IV 
ignore the effect of quark-pair creation (vacuum polarisa- 
tion) on meson masses. In this section possible modifica- 
tions to these predictions are considered in the framework 
of the unquenched quark model, which takes account of 
pair creation. 

In the strong coupling picture of a heavy quark QQ 
meson, the creation of a light qq pair is manifest as a cou- 
pling to a pair of heavy-light mesons Qq and qQ. Above 
threshold, this coupling leads to strong decay. Below 
threshold, it leads to a virtual Qq qQ component in the 
wavefunction and shifts the physical meson mass with 
respect to its bare mass. 

Most models for the coupling assume that the qq pair 
is created in spin-triplet: where models differ is in the 
treatment of the spatial degrees of freedom. In the con- 
text of decay this assumption first appeared due to the 
requirement that the created qq pair has the 3 Pq quan- 



tum numbers of the vacuum [57], but it is also follows 
from invariance arguments [58] and is a feature of the 
flux-tube model of strong-coupling QCD [59-61]. The 
first lattice QCD calculations of hybrid meson decay [62] 
arc also consistent with the spin-triplet hypothesis [63]. 

In refs [64, 65] a general expression is presented for the 
matrix element for the coupling QQ <-> Qq qQ in mod- 
els which assume spin-triplet pair creation. Barnes and 
Swanson [66] apply the expression in perturbation the- 
ory to compute the mass shift due to quark-pair creation. 
Their key result, which applies in a symmetry limit, is 
that the shift is independent of the spin and total angu- 
lar momenta of the QQ. Thus, for example, the masses 
of bare states n 1 So and n 3 Si are shifted equally and the 
hypcrfinc splitting is not affected. 

In practice this ideal limit is not realised: the sym- 
metry is spoiled by the different masses of the bare QQ 
states, and of the heavy-light states Qq and qQ to which 
they couple. Bare states QQ with different spin and total 
angular momenta are shifted by different amounts, lead- 
ing to induced spin-dependent splittings. The question 
relevant to this paper is how these effects modify the pre- 
dictions for IS and 2S hyperfine splittings. Refs [67, 68] 
address a closely related question, of how quark-pair cre- 
ation modifies the hyperfine splitting of P-wave mesons. 
It turns out the the same arguments carry over to S-wave 
hyperfine splitting, so the approach and its conclusions 
are briefly summarised here. 

In the non-relativistic limit of the quenched quark 
model (which ignores pair creation) the hyperfine split- 
ting 

i (M 3Po + 3M 3Pl + 5Msp 2 ) - Mi Pl , (20) 

is exactly zero, a prediction strikingly consistent with 
experiment for IP charmonia [69] and for IP and 2P 
bottomonia [70, 71]. The induced mass shifts of the un- 
quenched quark model might be expected to spoil the 
cancellation in the above linear combination, but remark- 
ably this turns out not to be the case. 

In refs [67, 68] it is shown that there is a mechanism 
in place which keeps the P-wave hyperfine splitting small 
despite large corrections to the meson masses. The proof 
considers the effect of coupling to combinations of heavy- 
light pseudoscalar and vector mesons 1 So and 3 Si in per- 
turbation theory. (The shifts due to coupling to orbitally 
excited mesons are smaller, being suppressed by an en- 
ergy denominator, and are in any case closer to the ideal 
limit of Barnes and Swanson [66] and so can largely be 
absorbed into a re-definition of model parameters.) As- 
suming only the creation of a qq in spin-triplet, the mass 
shift is expressed in terms of a power-series expansion in 
a parameter which is small provided that the mass differ- 
ences among the bare and heavy-light mesons are small 
compared to the binding energy. 

The result is that to first order in the expansion param- 
eter, the effect of quark-pair creation is to suppress the 
hyperfine splitting with respect to that of the bare states 



6 
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P Q<3 


AM4 S 


(am; s ) 


Barnes 


IS 


118 


0.69 


84 


(81) 


and Swanson [66] 


2S 


49 


0.51 


25 


(25) 


Kalashnikova [72] 


IS 


129 


0.899 


117 


(116) 


2S 


64 


0.743 


48 


(48) 


Eichten et al. [73] 


IS 


120.7 


0.966 


117.1 


(116.6) 




2S 


67.2 


0.702 


46.3 


(47.1) 



TABLE II. Hyperfine splittings AM nS and AM' nB of charmo- 
nia in quenched and unquenched quark models respectively. 
The final column (AM^ S ) is the approximate splitting ob- 
tained from equation (22). The probabilities Pqq in equa- 
tion (22) should strictly be computed using spin-averaged me- 
son masses, whereas those which appear in this table are the 
quoted probabilities for the n 3 Si states computed using phys- 
ical masses. 



by the factor Pqq, the probability that the physical me- 
son is in the bare state QQ rather than the meson-meson 
state Qq qQ 



QQ 



= 1 



P Qq qQ' 



(21) 



(Strictly, this is the probability computed using spin- 
averaged masses, though in practice it is similar to that 
using physical meson masses.) Consequently, if the 
bare states have zero hyperfine splitting as in the non- 
relativistic quark model, then so do the physical states, 
apart from small higher order corrections. In models in 
which the bares states have a non-zero hyperfine split- 
ting, such as those which include relativistic corrections, 
the effect of quark-pair creation is to reduce the splitting. 

For the hyperfine splitting of S-wave mesons, it turns 
out that the same result applies. The hyperfine split- 
ting AM^ S after including the effects of pair creation is 
suppressed compared to is bare value AM„s according to 



AM; S = P QQ AM nS . 



(22) 



The derivation, which is analogous to but more straight- 
forward than in the P-wave case, will be discussed in 
a future paper investigating generic features of the un- 
quenched quark model. For the purposes of the present 
paper, its validity is verified by testing it against a sev- 
eral model calculations in the literature. Three examples 
are shown in Table II. 

The first obvious consequence of equation (22) is that 
switching on quark-pair creation effects decreases the hy- 
perfine splitting, since Pqq < 1. This has immediate im- 
plications for the interpretation of the experimental data 
of Dobbs et al. and Belle. 

In Figure 1, all of the predictions for the 2S splitting in 
quenched quark model (with one exception) are smaller 
than that observed by Dobbs et al. Incorporating quark 
pair creation will suppress these predicted splittings fur- 
ther, exacerbating the disagreement. 

On the other hand, all of the predicted 2S splittings 
(with one exception) are consistent with or larger than 
that observed at Belle. In these cases, switching on pair 



creation has the potential to improve the overall agree- 
ment with Belle. 

Equation (22) also has implications for the ratio of 2S 
to IS splittings. Intuitively, one expects the probability 
Pqq to be smaller for 2S than for IS states, since the 
former are closer to threshold and therefore couple more 
strongly to heavy-light meson pairs. The examples in 
Table II are all consistent with this expectation. Equa- 
tion (22) then implies that switching on pair creation 
decreases the ratio of 2S to IS splittings, 



AM^g AM 2g 
AM{ S AMis" 



(23) 



As well as the examples in Table II, the inequality is 
satisfied by other model calculations [74, 75] which (be- 
cause the probabilities Pqq arc not quoted) have not 
been included in the table. (In the case of ref [75] the 
inequality is very close to an equality.) In the context of 
charmonia, refs [73, 76] note that the 2S splitting is sub- 
stantially reduced by pair-creation effects (compared to 
the IS splitting), and that this improves the agreement 
of their model with experimental data; this observation 
is consistent with the above. 

Again, the effect makes it more difficult to reconcile 
the observations of Dobbs et al. with theory. The ra- 
tios of 2S to IS splittings of all of the quenched quark 
model predictions collected in Table I and Figure 1 ex- 
ceed the observed ratio of Dobbs et al, as discussed ear- 
lier. Switching on the effects of pair creation will decrease 
the ratio further, worsening the disagreement. On the 
other hand, all of the ratios are consistent with, or larger 
than, the observations of Belle, so quark-pair creation 
could potentially improve the agreement. 

Finally, the effect of switching on pair creation on the 
on the model-independent relation of Section III will be 
considered. As discussed, the hyperfine splitting AM„s 
of a given level is suppressed with respect to the quenched 
quark model value by a factor Pqq- For the e + e~ width, 
the virtual photon couples to the QQ — as opposed to 
Qq qQ — component of the meson wavefunction. Ignor- 
ing the mixing between different n levels induced by pair 
creation, which is typically small in models (see for ex- 
ample refs [73, 77]), the e + e~ width of the unquenched 
quark model is related to that of the quenched quark 
model by 



r e+e-^n 3 Si ~ PqQ^c+c- ->n 3 Si • 



(24) 



The ratio (14) of hyperfine splitting to e + e width re- 
mains as it is in the quenched quark model 



AM, 



nS 



AM, 



nS 



r' + 

e + f 



>ra 3 Si 



r e+ 



e+e~— ¥n 3 Si 



(25) 



and the model-independent relation holds. 
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VI. CONCLUSION 



On the basis of their predictions for the hyperhne split- 
ting, each of the theoretical approaches discussed in this 
paper supports the interpretation of the Belle candidate 
as the r]b(2S). Among the lattice QCD results, there is 
only one quoted value which disagrees with the Belle hy- 
perhne splitting, and this is shown to be consistent when 
normalised to more recent, smaller experimental values 
of the IS splitting. In contrast, there is only one lattice 
result (with large errors) which is consistent with Dobbs 
et al, and this is also consistent with Belle. 



A model-independent relation, which is satisfied for 
charmonia and is expected to be more reliable for bot- 
tomonia, is in agreement with Belle but not with Dobbs 
et al. The relation is shown to be consistent with lattice 
results. 

Many potential models are consistent with Belle, and 
several others may plausibly be brought into agreement 
by adjusting a s or taking account of mass shifts due to 
quark-pair creation. On the other hand, models consis- 
tently predict the hyperfine splitting much smaller than 
that of the Dobbs et al. candidate, and it appears un- 
likely that adjusting model parameters or incorporating 
quark pair creation could bring them into agreement. 
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